Preparing high purity initial states for 
nuclear magnetic resonance quantum computing 
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Here we demonstrate how para-hydrogen can be used to prepare a two-spin system in an almost 
pure state which is suitable for implementing nuclear magnetic resonance (NMR) quantum compu- 
tation. A 12 ns laser pulse is used to initiate a chemical reaction involving pure para-hydrogen (the 
nuclear spin singlet of H2). The product, formed on the /is timescale, contains a hydrogen-derived 
two-spin system with an effective spin-state purity of 0.916. To achieve a comparable result by direct 
cooling would require an unmanageable (in the liquid state) temperature of 6.4 mK or an impractical 
magnetic field of 0.45 MT at room temperature. The resulting spin state has an entanglement of 
formation of 0.822 and cannot be described by local hidden variable models. 

PACS numbers: 03.67.Lx, 03.67.Mn, 82.56.-b 



Introduction. While quantum computing |l| offers the 
potential of using new quantum algorithms to tackle 
problems that are intractable for classical processors, 
its implementation requires the development of quan- 
tum devices, which are as yet unavailable. The most 
complex implementations of quantum algorithms to date 
have used techniques adapted from nuclear magnetic res- 
onance (NMR) spectroscopy 0,11110, but current liquid 
state NMR approaches cannot be extended to systems 
with many quantum bits, as it is not possible to prepare 
pure initial states by directly cooling the spin system into 
its ground state 0]. Furthermore, it has been shown that 
current NMR experiments involve only separable states 
, and thus could in principle be described by local hid- 
den variable models. 

The conventional approach in NMR quantum comput- 
ing is to use an ensemble of spins, and to prepare a 
pseudo-pure ground state d Q of the form 



Para-hydrogen induced polarization. A pure singlet 
nuclear spin state can be described using product op- 
erator notation [la as 
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where 1/2" is the maximally mixed state of an n-spin 
system, and e is the polarization of the state. In the high 
temperature regime this approach is exponentially ineffi- 
cient Hi . Furthermore if the polarization lies at or below 
a critical bound then any apparently entangled states 
prepared from the pseudo-pure state are in fact separa- 
ble 0; for two qubits [1 |l [lo| this bound is e = 1/3, 
corresponding to a fractional population of 1/2 in the 
ground state and 1/6 in each of the three other eigen- 
states. A radically different approach is to prepare initial 
states using non-thermal means [lll |. e.g., by using the 
pure singlet nuclear spin state isomer of H2 (called ^^para- 
hydrogen") |l2jllHlJ)ll3 ^ '^o^'^ spin-state reservoir. 



where E — I^ = ^(ca;®!), Sx = ^{ti^ax) and so 

on. The existence of para-hydrogen is a consequence of 
the Pauli principle 17], which requires the overall wave 
function of the molecule to be antisymmetric with re- 
spect to particle interchange. Dihydrogen molecules in 
even rotational states (J = 0, 2, . . .) possess an anti- 
symmetric nuclear wave function and correspond to nu- 
clear spin singlets (5*0 = (| Ti) — I iT))/V2 = termed 
para). Molecules in odd rotational states (J = 1, 3, . . .) 
have symmetric nuclear wave functions and consist of the 
three nuclear spin triplets (Tq = (| ti) + | iT))/^/2 = 
T-i = Itt), and T+i = |U), termed ortho). Note that 
the T±i states are not the same as = (| tT)±| ii))/'\/2 
(the other two Bell states), but that an equal mixture of 
T4-1 and r_i is also an equal mixture of and <I>~. 

Isolation of the para-hydrogen spin isomer is possible 
because spin-isomer interconversion is forbidden by angu- 
lar momentum selection rules. Adsorption onto a suitable 
surface breaks the symmetry of the II2 molecules, allow- 
ing spin-isomer interchange. The new ortho /para ratio 
therefore remembers the temperature of the last conver- 
sion surface encountered. Upon moving away from this 
surface interconversion is again suppressed. A tempera- 
ture of 20 K is sufficiently low to form the J = state 
and hence produce pure para-hydrogen. The resulting 
para-ll2 molecule cannot be used directly, as it is NMR 
silent due to its high symmetry. By means of a chemical 
reaction, producing a new molecule, the two hydrogen 
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atoms can be made distinct (I and S) and can be sepa- 
rately addressed. This phenomenon is well known from 
mechanistic NMR studies of catalytic hydrogenation and 
hydroformylation, where it is usually referred to as para- 
hydrogen induced polarization, or PHIP [ij, [ij, [ij, . 

In many previously described PHIP experiments the 
addition process is slow in comparison with the frequency 
difference between the I and S spins of the reaction prod- 
uct. This causes the off-diagonal terms in the density ma- 
trix to dephase, resulting in an equal mixture of singlet 
and To triplet states in the spin ensemble which is 
separable 0, 0, 0| . One possible method that has been 
suggested to overcome this problem is to perform the slow 
addition while applying an isotro pic mixing sequence [Tflj ] 
to remove the dephasing effect |2(J|. In principle (and 
neglecting relaxation effects) this should completely con- 
serve the hydrogen spin state, and this method has been 
used [23 to achieve states with e « 0.1. 

A much simpler approach is to ensure that addition is 
rapid in comparison with the dephasing and relaxation 
timescales. This requires addition to a highly reactive 
species, but a fast reaction would go to completion and 
the spin system would decohere before the sample could 
be studied. In this paper we solve this problem by us- 
ing photochemistry to prepare a reagent which will re- 
act instantaneously with paro-hydrogen to produce the 
molecule for examination. This means we can start and 
stop the formation of the reactive species in a controlled 
manner. We then show that the resulting two-spin sys- 
tem forms in an almost pure state. 

Experimental methods. In this study the precursor 
molecule was Ru(C0)3(dppe), where dppe indicates 1,2- 
bis(diphenylphosphino)ethane. This was prepared from 
Ru3(CO)i2 by warming a benzene solution to 373 K un- 
der 30 atmospheres of CO in the presence of three equiv- 
alents of dppe 21]. The Ru(C0)3(dppe) was dissolved in 
dg-benzene in a 5 mm NMR tube fitted with a Young's 
valve to permit attachment to a vacuum line. Dissolved 
gases were removed by freeze-pump-thaw cycles and the 
tube was covered in foil to exclude light. Para-hydrogen 
was prepared at a temperature of 20 K using a charcoal- 
based ortho-para interconversion catalyst; at this tem- 
perature the thermal state is essentially pure para. The 
para-hydrogen was then introduced to the sample; as the 
catalyst is no longer present ortho-para interconversion is 
suppressed. After warming, shaking ensures that the II2 
gas (with a pressure of about 3 atmospheres) dissolves. 

The NMR tube was placed in a 400 MHz NMR 
spectrometer fitted with a ^H/^^P tuned NMR probe 
equipped for in situ photolysis |2^ . (The transfer was 
performed in a darkened room to prevent premature pho- 
tolysis by ambient light.) The reaction was initiated by a 
12 ns pulse of 308 nm UV light (pulse energy 32 mJ) from 
an MPB Technologies MSX-250 pulsed XeCl excimer 
laser triggered by the NMR spectrometer. This generates 
the reactive intermediate Ru(C0)2(dppe), in situ from 



its stable precursor Ru(C0)3(dppe) by laser flash photol- 
ysis (see Fig. 1). Subsequent reaction of Ru(CO) ^(d ppe) 
with H2 occurs on the sub-microsecond timescale j23] and 
leads to the product of interest, Ru(H)2(CO)2(dppe). 

The two hydride resonances of the product appear at 
S — — 7.55ppm (spin I) and (5 — — 6.32ppm (spin S). For 
the analysis a spin I selective pulse was implemented us- 
ing two hard 90° pulses separated by a delay of 1/(45:^), 
where — 492 Hz is the difference between the reso- 
nance frequencies, with relative phases of 135° and with 
the RF frequency centered on the midpoint of the two 
resonances. Such pulses, based on Jump and Return se- 
quences, have been described previously ji^l. The GARP 
sequence |2^ was applied throughout signal acquisition 
to remove couplings to '^^P nuclei. To avoid the necessity 
for complete quantum state tomography |26l | a filtration 
sequence was developed which has no effect on the de- 
sired singlet state, but dephases most other states. 

While the signal produced by the pora-hydrogen is eas- 
ily seen in one scan, the thermal signal is extremely weak 
and so difficult to measure directly. It was therefore nec- 
essary to increase its intensity by increasing the amount 
of Ru(H)2 (00)2 (dppe) in the sample. This was achieved 
by applying a further 999 laser pulses to produce more 
Ru(H)2 (00)2 (dppe). Even then we needed 3072 scans 
for the calibration spectrum; these were separated by an 
interval of 20 s, which is much greater than five times 
the measured Ti of 1.7 s and so saturation effects can 
be ignored. The spectra were processed by homewritten 
software and analysed by integration. The para-hydrogen 
spectrum comprises a pair of antiphase doublets, which 
partially cancel jl^], and so direct integration will result 
in an underestimate of the signal intensity; to reduce such 
effects the spectrum was J-doublcd 28] four times before 
integration. Even after J-doubling and integration slight 
imbalances were visible between the two multiplets, and 
these imbalances can be analysed to determine the im- 
balance between Tq and T±\ triplet states. 

Results. Laser photolysis of the Ru(C0)3(dppe) leads 
to the product Ru(H)2(CO)2(dppe). If the two hydrogen 
nuclei have inherited the nuclear singlet from the para- 
hydrogen, then a selective 90° ly pulse will yield the ob- 
servable NMR terms 

i(-2/,5, + 2/,5,). (3) 

This corresponds to a pair of antiphase doublets, with 
intensities ±1/2. This pattern is indeed seen (Fig. 2), but 
to show that we have an essentially pure singlet state it is 
necessary to determine the intensity of the signal as well 
as its form. We calibrated our signal against a standard 
provided by the thermal state of the same spin system, 
obtained by allowing the spin system to relax. At high 
temperatures the thermal state is 

\(\E^\B\h^S.\) (4) 
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FIG. 1: The reaction scheme used to generate Ru(H)2(CO)2(dppe) with an almost pure initial spin state. A UV photon knocks 
one carbonyl group off the Ru(C0)3(dppe) precursor to generate an unstable intermediate, which immediately adds hydrogen 
to give the desired product. Since para-Yi.2 has a pure singlet initial state, and addition occurs with retention of spin state, the 
product should also have a pure singlet spin state. 
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FIG. 2: The para-hydrogen enhanced and calibration spectra 
of Ru(H)2(CO)2(dppe). The para-hydrogen spectrum (a) is a 
single scan after a single laser flash, while the calibration spec- 
trum (b) is the sum of 3072 scans after 1000 laser flashes. The 
calibration spectrum has been divided by 3072 x 1000 and then 
multiplied by the theoretical maximum enhancement (31028, 
see main text) before plotting, so that the two spectra should 
show the same intensity. In fact the para-hydrogen spectrum 
(a) is even more intense than naively predicted. 



where B — hv /kT and v is the Larmor frequency of the 
spins. After a 90° pulse the NMR observable terms are 



\B{h + S,), 



(5) 



that is a pair of inphase doublets with intensity B/A, 
so that a pure state will give a signal 2/B times larger 
than a thermal state. For our system v = 400 MHz and 
T = 295 K, and so the para-hydrogen signals could be up 
to 31028 times more intense than the thermal signals. 

In order to examine the form of the signal, we used a fil- 
tration sequence (comprising two periods of length l/Sv, 
under a field gradient separated by a hard 90° pulse) 
which is closely related to "twirl" operations [23. l30| . 
After this, the density matrix is Bell-diagonal, compris- 
ing a mixture of singlet and triplet states, with equal 
amounts of the parallel-spin triplet states T+i and T_i. 
Since this sequence had little effect on the spectrum, the 
initial density matrix had a similar form. 

Our results (Fig. 2) show an apparent enhancement 
of about 77000, significantly higher than expected. This 
discrepancy arises because the NMR probe is not sensi- 
tive to the entire sample, but only to that within the RF 



coil. The hydride forms within the coil region because of 
the position of the mirror that introduces the UV light, 
as confirmed by one-dimensional NMR imaging, but is 
then distributed throughout the sample by convection 
and diffusion, so that in the calibration spectrum only a 
fraction of the hydride is detectable. As expected the di- 
rectly measured enhancement (that is, before correction) 
shows a linear dependence on the total sample volume 
(data not shown). The active volume fraction can be esti- 
mated using geometrical arguments based on the relative 
length of the NMR sample and the RF coil. The actual 
value depends on the length of the NMR sample, but was 
0.368 in the experiment described. After correcting for 
this active volume fraction the observed enhancement is 
consistent with a polarization of e = 0.916 ± 0.019. 

Relaxation. The measured relaxation and decoher- 
ence times of the nuclei in the para-hydrogen state 
(Ti = 1.7 s, T2 = 0.58 s) are indistinguishable from those 
in the thermal state, indicating that the high polarisation 
does not affect the relaxation properties of the molecule. 
No effects of radiation damping were observed, reflecting 
the extremely low concentration of the hydride. 

Entanglement. As well as the polarization enhance- 
ment, the use of PHIP initializes the system directly into 
an entangled state. The entanglement threshold of 1/3 is 
valid only for mixtures of a Bell state with the maximally 
mixed state (Werner states). Our density matrix could 
contain arbitrary states mixed with the singlet and the 
threshold depends on what is mixed in. The states that 
most effectively destroy the entanglement of the singlet 
are Tq and equal mixtures of T-|_i and T_i; these are also 
the only states which survive the filtration sequence. 

The positivity of the partial transpose test |3| is a nec- 
essary and sufficient condition for separability for two 
qubit states P|. It can be shown that for a sys- 
tem comprising a mixture of the singlet and some convex 
combination of Tq and an equal mixture of T+i and T_i 
the state will always be entangled if the total amount of 
singlet exceeds 1/2. Detailed analysis of our data shows 
that the density matrix is not quite a Werner state: in- 
stead the state comprises 93.7% ^o, 4.5% Tq, 0.9% T+i 
and 0.9% T-i. The excess population of Tq suggests the 



4 



presence of some phase decoherence process that we have 
yet to identify, but this excess has no effect on the en- 
tanglement of formation [2^ l3]| , which would be 1 for a 
pure singlet and for our system is 0.822 ± 0.039. 

Our results demonstrate that liquid-phase NMR can 
produce entanglement. The question of whether entan- 
glement really is a necessary resource for universal quan- 
tum computers can only be answered by mathematical 
proofs of their scaling characteristics with and without 
entanglement. These can only be defined rigorously in 
the limit as the size of the problem instance (and hence 
the computer solving it) tends to infinity. Any actual 
experiment can only be performed on a finite-sized sys- 
tem, and can therefore neither prove nor disprove the 
correctness of the arguments in [7|. We can only show 
that this particular objection can no longer be levelled at 
liquid-phase NMR quantum computation. 

Conclusions and further work. We have shown that 
para-hydrogen can greatly benefit NMR quantum com- 
puting by the generation of almost pure (e = 0.916 ± 
0.019) initial states on demand, without the need for 
lengthy preparation sequences. To achieve a comparable 
result by direct cooling would require an unmanageable 
(in the liquid state) temperature of 6.4 mK or an imprac- 
tical magnetic field of 0.45 MT at room temperature. 

Conventional liquid-phase NMR has poor scaling char- 
acteristics at low polarization 1^, but our effective spin 
temperature of less than 6.4 mK is much less than the 1 K 
threshold suggested by Warren for efficient initialization. 
This approach can in principle be scaled up by adding M 
molecules of para-hydrogen to a single precursor, effec- 
tively synthesizing a quantum computer with 2M qubits 
in a pure initial state, and we are currently seeking to im- 
plement these ideas. We also expect that the ability to 
rapidly generate highly spin-polarised hydride species on 
demand will prove useful in other para-hydrogen based 
studies of reaction mechanisms. 
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